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SUMMARY 


Ar. analytical model consisting of two ring sources Is developed to 
study the direct radiation of jet noise In te ms of correlation, coherence, 
and phase and also to aid In solving the inverse radiation problem of 
determining the no<se source In terms of far-field measurements. The rings 
consist of discrete sources which are either monopoles or quadrupoles with 
Gaussian profiles. Only adjacent sources, both within the rings and between 
rings, are correlated. Results show that from the far-field Information one 
can determine when the sources are compact or noncompact with respect to 
the acoustic wavelength and distinguish between the types of sources. In 
addition, from the inverse radiation approach one can recover the center of 
mass, the location and separation distance of the ring and the diameters. 

INTRODUCTION 

This paper addresses a solution of the direct and Inverse problems in 
the radiation of sound from randomly fluctuating sources on two rings. This 
effort resulted from the fact that even a well designed experiment cannot 
be free of ambiguity and as such it Is shown that a parallel study using a 
model Is essential in minimizing the experimental uncertainties. In addition, 
the recovery of the unknown source distribution on a jet is more easily 
obtainable from an analytical model, which is by far a superior tool th^n 
from the measurements. 

The approach utilizes space-tiiie correlation of the far-field sound which 
clearly provides more information of physical importance than a single point 
data. The first theory and experimental work on two point correlation from a 
jet were reported in references 1 to 3, and for the inverse approach in 
references 4 and 5. Several other papers have been written discussing the 
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recovery of the unknown source field In an application to the jet noise 
problem. The most recent ones are given in references 6 to 10. 

Conceivably, there are n.any possible ways to utilize the radiation 
data to yield some Information about the source field. In general, however, 
the common difficulty In dealing with Inverse problems In physical sciences 
lies In the fact that they are Ill-posed In Hadamark sense (ref. 11). The 
present approach was designed to fit the physical problem under consideration; 
that is, treating the direct and Inverse stochastic problem In a uniform 
media for jet application. 

The first part of the analysis contains the direct approach and numerical 
results on two ring sources. The sources In the rings are random and coupled 
with Gaussian profiles. To better distinguish the type of singularities, 
two types were used; (a) random monopole and (b) random quadrupole. The 
second part of the paper contains the inverse approach accompanied by the 
numerical results and the description of the method to determine the basic 
geometry with compactness and identification. 

ANALYSIS 

The source model consists of two rings, centered on and, perpendicular 
to the jet axis. Each ring is made up of a finite number of discrete sources 
which may be either monopoles or quadrupoles. The rings have different radii, 
simulating the spreading of the jet. The sources in the rings are random in 
time with Gaussian profile. Only adjacent sources, both within the rings as 
well as between rings, are correlated. 

In order to investigate the phase effects arising from different types of 
sources, two types of sources are used: (a) random monopoles and (b) random 

quadrupoles. 
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The variables In the model are the number of sources In the rings, 
distance between rlnqs, and rlnq radius. The sources can either be real or 
complex. A simple Interpretation for complex sources In the direction of 
the jet can be given In terms of convection effects. The sources, however, 
can also be Interpreted as the Llghthlll sources, that Is, the sources of an 
equivalent medium at rest. 


A. Direct Radiation 

The geometry of the model is given In figure 1. The symbol j(t) 
denotes the source placed at an angle 4>j In the 1-th ring with 
<t>j * (j - 1 )A, j * 1 . . . N and A * 2n/N. N Is the number of sources 
which are equally spaced over the rings of radii a^ for 1 ■ 1, 2. The 
ring with radius a ? Is larger than the ring with radius a^. The two 
rings are placed at distance and respectively from the exit of an 
ideal jet. 


Position 

Strength 

x y 

z 


1st Ring Sj jUj, aj cos j ♦ 

a 1 sin $■) 


2nd Ring Sp ^2 0 j » 

a 2 sin 4>j) 



£ k denotes the far-field position vector of the field points with k = 1, 2. 
In spherical coordinates, C^(r^, , ♦•j ) and Lpfr^, e^, We set 

r l * r 2 = r ’ case d ^ stance between a point on the ring and a 

point in the far-field Is given by 
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l**k ’ s ljl " y** 2 ♦ ♦ a? - 2rl ^ cos 9^ - 2ra^ sin 0^ cos(<>j - $ k ) 

( 2 ) 

- r - ^ cos e k - a^ sin 0^ cos (<J> j - $ k ) 

We denote the self correlation of the random source by (x ) and Q 2 (t). 
Therefore, 

p 2 < q l,j^ q l,J^ + T ) > * Q]( T ) 

P 2 < ( ^2,j^ t ^ q 2,j^ t + T ) “ ^2^ t ^ 


We assune that only the adjacent sources are correlated: 

p2 ' q i.j< t)q i.k (t + t) » ■ Q i.i (T) [ s i,u-k|] 

p2 * q 2„| ( ‘ )q 2,k (t * T > > ■ Q 2,2 (T >[ 5 MJ-k|] 

P 2 < q l f j^ q 2,J^ t + > s Q 1,2^ t ^ 

P 2 < q l,j^ * T ^ q 2,j^ > s ^l ,2^” t ^ 

p2 * q l,J (t)q 2,k (t + T > » ’ Q 1.2 (T) [ S l,|J-k|] 

p2 < q 2,k |t)q l,j (t 1 t) > ’ ’ , *,2 < " T) [‘'l,|j-k|] 

We implied that q^ q = q. and q^ ^ = q^. Note that 
Q- , are even functions of t. 


Q 1 ’ Q 2’ Q 1 ,1 and 


Is 
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The 


pressure P^ due to source ^ 


P. - pq 


(‘ - 


-)/■ 


|L - s < A 


where q^ j(t) ■ q^ j. The time derivative Is used to convert Into 
unit of pressure. The far-fleld cross correlation function Is: 


<P 1 (t)P 2 (t * t) > = p 7 < 


2 N 


h jl, "'-J 


{ . * 4 ' 4 .J 1 ' 


,|L t - s u 


2 

l 

k= 


l m (t 

1 m=l k,m \ 


t + T - 


^4 ' s k,nJ 


'14 ' S k.ml 


2 


2 N c. in / 

l l ^ j(t) l l q k m ft 

M j-1 1,J k-1 m=l k,m ' 

"•" re “l.J.k.n,* |-l l l • S 1.jl * 14 - S k,ml|/ c o 

cos 6^ + a f sin 0 ^ cos(0j - Jc Q 
-|[t k cos 3 2 * a k sin d., cos(* m - * 2 )]|/c 0 


2 N 


+ T - U 


1 » j * k ,m) 


> (3) 


(4) 


1. Rings with monopoles .- The far-field cross correlation for monopole 
sources can be written as follows 


< P,(t)P,(t ♦ T) > - - 4-4r, ? (t) 

f dt '•* 


^ \ w 2 
P 


( 5 ) 
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where 


( T • “lj.lj) * 

N 

A 

Q l.l( T ' W 1J,1 J+l) + Q l.l( T * w lJ,1.J-l) 

/ \ 

N 

* / \ * / \ 

( T ’ u lJ.2j) + 

l 

j«1 

Q 1,2( t ‘ w lj,2,J+l) + Q 1 .2 ( T ' u lj,2,j-l) 

/ \ 


Nf./ \ A / 

,2(‘ t + U 2J.U > 

1 + j«i p . ? V T ’ • 1 *J + V + Ql * 2 v T + ^ * ] *J 


.{ } Q 2( t ‘ p 2j,2j) + [ Q 2,2 ( T ‘ p 2j,2,j*l) + Q 2,2( T ' w 2j,2,j-l) 


We make use of the relationship 

d _ h t .\-/ \ ^ _ d 


< q(t)q(t + t) > * jp < q(t)q(t + t) > = - < q(t)q(t + t) > 

2 2 
* < q(t)q(t + t) > * - Q(t) 


( 6 ) 


dx 


d T ‘ 


The Fourier transform of ? {r) is 

A a N IdJUi i •) t 

R 1 > 2 (w) - Qj ('u)J e 

A N i U)U A J A J 

♦ Q 2 (u)J e 2 ‘ J ’ Z ’ j 


Uncorrelated sources 
on the first ring 


Uncorrelated sources 
on the second ring 


Q] 

1 » 1 i 


e 1 < 4 J l,j,l,j+l + e 1 wV 1 l . j +1 J . j 


Correlation between 
S l,j ; S l,j+1 




♦ Q?, 2 (“)j 




a N 1 U^J 14 0 4 A 

* Ql ,*<->{ e + Q 1 .2 ( ” u,) f 


N IlilVl 


2.J.1.J 


Correlation between 
S 2 , j ’ S 2,j+1 


Correlation between 


N Iwvi 


* lUHi| j f) j.* a. il 

0j i2 (-) { e 1Jt2 * J+1 ♦ 2 (-u) I 


N lu'ii 


2.J+1.1.J 


S U ; S 2,j 


Correlation between 
S 1J ; S 2,j+1 


N Iwu 


N lu'ii. 


♦ 6; 2 W I e * q; >t M i « J Corre'iilon 

' ’ S I,J : S 2.J-1 


between (7) 


2 2 2 ' 

For the computation of the cross spectral density, -(r /p ui )Rj ? (ui), one 
needs the six correlation functions of the sources; l.e., Q 1 , Q^, Q 1 y 
Q 2 2 * Q] 2’ and ^1 2 or the ^ r Fourier transforms. 

For tne purpose cf computlnq the cross correlation R^ 2 ( r ) let us 

assume that (^(t), Q 2 (t), Qj 1 ( t ) , ^2,2^ ^ * ^1,2^’ and ^1,2^ T ^ have 
Gaussian form which has frequently oeen used in the past (refs. 12 and 13). 


Qj (t ) “ e 


-( Yi t) 4 


. / 




Q 2 (t) * e ^ 


Ql,l (T) ’ 2 


U ■ , \2 

2 / TT 
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\2 

‘M T ‘ fi *i) 
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(t) 


Q 1.2 (t) 



(t) 





( 8 ) 


The constants y's control the rates of decay while i^'s represent the 
time shifts. The Fourier transforms of the Q's are then given by 



j(us) * Q^'uj) cos <1^) 


^2,2^^ * $ 2 ^ cos (lT ^ 2 ) 



( 9 ) 
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In accordance with our assumption of axl symmetry, we set "0*, 
and chose the fixed point to be at ■ 90*. The computation of Rj 
In the meridian plane with * 0* and In the azimuthal plane with e ? ■ 90* 
are shown In figures 2 and 3. The number of sources chosen In each ring 
Is 16. T he location of the rings are at ■ 2D and ■ 5D, D Is the 
diameter of the Ideal jet taken ■ 0.06 k ; 5 m corresponding to the diameter 
of the experimental jet. In addition, the spreading angle connecting the 
two rings of radius a 1 ■ .8460 and a ? ■ 1.36D was chosen to be 10° 
corresponding to a jet spreading half angle. 

The broadband cross correlation function (fig. 2) shows the decay 
and the shift of the peak correlation away from t ■ 0 with the increase 
In separation angle from 02 from 0^ = 90°. The same Is observed in the 
experimental data (ref. 2). The azimuthal broadband cross correlation 
function, figure 3, shows two results at ♦« ■ 180° corresponding to two 
different decay rates. One Is for a slow decay rate; the sources are compact 
as Is shown by the fact that the cross correlation peak at zero time delay. 
This is also true for $2 < 180 * The second result Is for fast decay 
rates; the peak cross correlation occurs at a time delay different from 
zero. Indicating that the sources are not compact with respect to the 
acoustic wavelength. In this example, the noncompactness contribution 
became noticeable only at large separation angle as thp cross correlation 



./ 

! 

# 



10 

approaches 180* separation In a broadband sense. Also, In th*. experimental 
results* the peak correlation occurs at x t 0 at large separation angle 
$2 - showing the nonco.-ipactness of the sources. The fact that the 
broadband cross correlation Is even In t Indicates that statistically 
the sources are symmetrical ly placed with respect to the axis of the jet. 

The phase of the far-fleld pressure In the meridian plane Is shown In 
figure 4. Results are plotted In terms of the frequency w for 
ranging from 80* to 0*. Note that the phase goes to zero with the frequency 
for all angular separations. This proves that the sources are monopoles, 
since no net cancellation occurs at any one angle more than at any other as 
u> goes to zero. Also, note that the Initial phase slopes are constant at 
low <*> but varying at higher u>. This u an Indication that at lower fre- 
quency the radiating field pattern Is nearly uniform, while at higher fre- 
quencies It Is not because the sources are noncompact. From this figure, 
however, one cannot deduce a priori that the sources are noncompact beyond a 
certain range of frequencies, without first establishing the phase due to 
geometrical effects. This problem will be discussed further as part of the 
inverse problem In section B. In addition, the pattern of the phase varia- 
tion with id and 0 resembles the measured phase (ref. 2), but with one 
Important difference; namely, that in the experiment the phase for less 
than 50° Is nonzero as the frequency approaches zero. This Is evidence that 
the sources In the experimental jet are not monopole. This point also will 
develop further when comparison is made between the phas. of monopoles and 
quadrupoles In section A-2. 

‘Results not yet published show that by filterinq out the lower frequency 
on the original broadband spectra, the peak correlations shift near'.y 
symetrlcally away from t * 0. 
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2. Rings with quad rupolef The far-fleld quadrupole of the pressure P 
In terms of the far-fleld anqle 0 and # and distance r are obtained as 
linear combinations of the fol lowing 
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P - 




2 

COS 0 


rfxdy 


P . 
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P - 


dxdz 


dydz 
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P . 


P . 


ccs 0 sin 0 cos 4> 


2 2 
sin e cos ♦ 


sin 0 cos 0 sin <l> 


sin 0 cos $ sin $ 


2 2 
sin 0 cos $ 


( 10 ) 


where x * r cos 9; y • r sin 0 sin <*; z ■ r sin 0 sin 1> 
The total contribution of a quadrupole becomes: 




2 2 2 
a^ cob 0 + aj 0 cos 0 sin 0 cob J> + a.,., sin 0 cob $ 


2 2 2 

+ a 13 cos 0 sin 0 sin 4> + a 23 sin 0 cos * sin ♦ + sin 0 sin <t 
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• - ^*4 (t " T" ) P(®» ♦) (11) 

c^r dt J \ c o / 

We will assume that the coefficients a^ , a^ a,, are 

Independent of t for all quadrupoles on the r<ngs. 

We make use of the following relationship: 


2 


P 




(t ♦ t) > 



( 12 ) 


The coefficient a 1 . . . . a 33 In F(G, $) will be adjusted to match 
with experimental results. The general form of the quadrupole correlation 
function, equivalent to equation 5, for the monopole becomes: 


r* 

~Z 

P 


< p,(t)p ? : 


. ) > 


r 2 

7 


d"7 R ’.2 (T) 


F(e r O 1 )F(0 2 , ^ 2 ) (13) 


where R^(t) Is defined In equation 5. 

Computation from the quadrupole model of the phase a and 

A 

| 2 ( u, )I/ r i ] (u) using the same meridian plane of the previous model Is 
shown '•> figures 5 and 6. The phase of the quadrupole differs from that 
of the monopole by 180* for * 45* as the frequency w goes to zero. 
Therefore, the type of singularities radiating sound in the far-field are 
distinguishable by the behavior of a as w tends to zero. At higher fre- 
quencies one cannot. In general, recognize the type of singularities because 
of their mutual Interaction. 
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The coherence function for the quadrupole model Is shown In figure 6 tor 
four different frequencies. It Is clear that at lower frequencies the pattern 
of the quadrupoles similar to the measurements made for an a xy symmetric 
Jet In reference 3. Indeed the phase toqether with the coherence function 
shows what types of quadrupole (namely the coefficients a^ . . . . a^) 
dominate. As the frequency Increases the coherence function decays due to 
randomness. 

The data shown In the figure have been obtained when all the correlations 
function Qj , . . . . Q^» are even In t. It is known, however, 

that, at least In the plane of the jet along the direction of flow the source 
cross correlation function Q^, Qjp are not even because of convection 
effects. Tills additional effect will not be discussed at this time. The 
present paper attempts to Interpret the results in this simplest form, that 
Is, to study the phase variation due to geometrical effects and compactness 
for different types of singularities. 

In order to study further the relationship between the quadrupole model 
and the measurements, comparison Is also made of the azimuthal cross power 
spectral density. Only the real part Is used since in the model the Imaginary 
part is very small at relatively lower frequencies. The calculated real part 
of the cross spectral density taken at 0^ * 0, * 90°, and ^ • 0° is shown 
In figure 7. The behavior is consistent with the measurements. Indicating 
the quadrupole model may be satisfactorily used to model the singularities 
in a jet. 
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B. The Inverse Radiation Problem 

A question naturally arises regarding the possibility 01 uslnq acoustic 
radiation data to infer the characteristics of the source distribution. 
Although the difficulty In dealing with this, and other types of Inverse 
problems In physical sciences, lies In the fact that they are In general Ill- 
posed, there may conceivably be some way to utilize the singularity model 
developed together with the measured radiation field to study tie Inverse 
relationship. 

Since there are many parameters In the model solution, we examine the 
solution In the meridian plane, that Is 4^ * 4\, * 0° and 0^ ■ 90*. 

Equation 7 reduces to: 


R 12 (w) - e 


-Ik^cos 


a H -1 cos 4 > . (A, ♦ B, 

Q, M r e j 1 1 

1 J-l 


+ Q,(u>) l 

' J-l 


-1(A 1 cos - B 1 cos 


-1 (A 1 cos - B 1 cos .j ) 


a N -1 cos <j> .(A, - B, ) 

♦ Q, 2 W l e J 2 1 

j*l 

a,. N -1 (A 0 cos 4 ) 4 ., - B, cos <j>.) 

♦ q: 2 m l e 2 j+1 1 j 

x%i J »1 

\[^2 cos <fj_j - B 1 cos 
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-Ik 1, cos 0, a N -1 cos (A 9 - B,) 
♦ e 0 2 2 Q>) l e j 2 ? 

2 J-l 


♦ Q ? ? («*») I e 
Z,Z j“l I 


N -1(A 2 cos <t>j - Bp cos $j + ]) 


-1 (Ap cos <J>j - B ? cos <f> j -j ) j 


where 


N -1 cos 4 j ( A 1 - Bp) 


♦ Qi o(w) l e 
},£ j«l 


a N -1(A, cos - B- cos <t > 4 , 

♦ Q? 2 Cu.) I e 1 32 J -' 

U£ J-l 


-1(A ] cos - Bp cos (14) 


k 0 a r* A 1 * - a 1 k o’ B i " - a 1 k o sin °2 ; 1 B 1 * 2 
0 

The second step consists of trying to simplify the above equation by observing 
the behavior In three domains: 1 ) k Q a i « 1 , 11 ) k o d^ 1 , and 1 i 1 ) k Q a i » 1 

1. Geometrical ce nter of t he sin gularity k Q a. « 1. - At large wavelength 
1i the limit as w -*■ 0, A^, B^ both are « 1, the first sur*. in equation 14 
reduces to 


t 

1 
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-Ik l. cos 0, a 
e 01 2 Q,(u>) 


N » 1 (A, - B,)^ COS *j * 0 [( ko a,) ? ] 


-Ik COS 0, * 

■e 0 ' 2 <),(«,) 


N * 0(k o a, )' 


II 

because ) cos • 

J»1 J 


Using a similar approacn for the remaining sums In equation 14 and omitting 
terms of 0 (( *c Q a ^ ) 2 ) we have 


R 1f2 M - N 


L -Ik 
(q^uje 


1., cos 0 O a -Ik?- cos 0 O 

0 ' 2 ♦ q 2 ( U )e 0 2 2 


1 


(15) 


where 


Q^(w) * Qj(w) + 2Q 1 j(w) + Q^ 2 (-u>) + 2 (-o)) 


a » % a ^ ^ 

q 2 ( u >) - Q 2 ( w ) + 2^2 2^ + ^ 1 , 2 ^ + *^ 1 , 2 ^ 


(16) 


If we also assume k Q fj « 1 and neglect the second order term we have 


~~2 2 («) “ n [R] (“)(! ” cos 0 2 ) q 2 (i»>)(l “ 


1k„fc 0 COS 0 O 
0 Z Z 


>] 


- N^q-j (uj) + q 2 (w)Je 


-Ik cos 0 O A 
0 z 


2 A 

~2 I Ri 2 vw) I e 
P C '•* 


laid) 


( 17 ) 




^iQi(w) + i * r* » 1 

A ■ — , a, ■ cos 6-A, |R, ~(w)| ■ N q, (u) ♦ q ? (u>) 

q^u.) ♦ q 2 (u>) 1 c o 2 1.? LI ? J 


The far-fleld phase Is Independent of u> correspondlnq to the geometrical 
center of the source (or center of mass). It Is the- Dhase due to the geo- 
metrical position between sources and far-fleld measurements. When A • 0, 

Is zero. 

A 

The magnitude |R 1 2^)1 as we H as A are available from measurements. 
The behavior of A Is shown In the sketch. The minimum value corresponds to 
the fixed microphone located at 0^, the center of the sources. 


Sketch - Behavior of the 
geometrical center for fixed 
far-field angle differences. 
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In fact, the results show that one can eliminate the phase resulting from 
the qeometry of the orientations. A numerical example Is shown In figure 8 
where the phases due to monopoles are evaluated with respect to the geometrical 
center. When the phases are evaluated differently, for example, with respect 
to a fixed point in the plane of the Jet (0^ * 90* shown In fig. 4), the 
results are radically different than that of figure 8. From the geometrical 
center the phase at low frequency Is zero as predicted by equation 17. The 
remaining phases as the frequency Increases are due to noncompactness of the 
source, l.e., ring separation, ring diameter, and due to the spacing between 
sources on the ring. 

Equation 17 describes the geometrical center only and cannot provide 
additional Information than what has already been discussed unless 
n,(u>) ■ q 2 (w); In general, however, this Is not the case. The geometrical- 
center result per se Is Important, especial ly from the measurement point of 
view, because the far-fleid distance can be taken from this center rather 
than using other center references. In so doing, the interpretation of the 
results will be considerably simplified. To obtain additional information 

regarding the source field, from far-field measurements, one can expand the 

2 

exponentially function by retaining terms in w rather than In w only. 

Such an approximation will be valid for k Q a. - 1 and is carried out in the 
next subsection. 

2. Ring location for k p a^ - 1.- Let us expand equation 14 In w and 

2 2 2 

retain terms of order u> or k Q a^ . The first term in equation 14 contains 
Q 1 (to) , can be written 
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Qi(w) I 
J-l 


\ cos (u>) - A ? (ui)j 


where 


Qt ( a>)N - n } k c 


Q*l (w)Nd« P 

n, « — - (1 - sin 0 2 r 


N N 

since [ cos 2 * 1 [ (1 ♦ cos 2 ■ 5[ 

j-1 J 2 j«1 J 2 


Following the same procedure for the additional seven terms In equation 14, 
one obtains the following equations for the cross spectral density: 


2 „ -Ik l. cos 0, - 2 

^R 12 (u.) . qi (u,)e 01 2 0 


+ q 2 (u)e 


-1 k L cos 0, - fl^k 2 

O 2 2 0 


(18) 


where 

* f2-| (ui) ♦ 2n 1 i (co) + 2 (-u») + 251* 2 (-a>) 

^ * n 2 (u)) + 2n 2j2 (o)) + n 1>2 (uj) + 2n * gM 

Qi ^ (w)Na£ « 

llj ^ L (1 *■ sin 0 2 •• 2 cos A sin © 2 ) 

Qi o(-w)N ? 

n i ,2^ _u> ^ " — 4 ^ a 2 ■ a 1 sin e 2^ 
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Q* ^2 + a 2 $1n 2 ^ . 2 cos A a,a 2 s1 " ‘ ^ 


T 

.2 


Q.(ii))N3a . »2 

fl 2 <»> ■ J T - <’ • 5l " e 2 ’ 


n 2 , 2 (u.) 


Q- ™(w)® 2 ^ , 2 ft \< 

. 2 » z . — t- (1 ♦ sin 0 2 ) 

4 


and 


Q^UON 2 ( 2 s(n 2 - 2 cos A V ? 5 '" 6 ' 
fl 1 2 (“>) ■ ‘"4 '“1 2 c 


Q. «(u>)N \! 

n, 2 W • M — (a 2 ‘ *' 2 


Equation 18 can be written as 

■Vi cos 9 2 [r;°)m^(vi^! 2)(u) ] 

) 2 R^ (<*>)] 


*7 * 1 . 2 ^ " 6 

P 


* e -V2 005 02 [i‘°>(.) - i (V2> 2R 2 


(19) 


rJ°'( U ) ■ q, Ml 4 ° >( “> ’ 

i‘ 2 >w ■[«,(») * * 

• COS T 5 ,ne 2 ) 


sin P 2 V 


where 


rj 
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♦ [q 1#2 (-w) ♦ 2 5l,2 t " w) ](a“ * s1n °tj 2 

r| 2) (u») • [q 2 (o>) ♦ 2Q 2t2 (a.)J(l - sin e 2 ) 2 

+ [^2,2 (tl,) + a" 2 ' cos t) s1n 0 2 

+ [^1.2 (w) * 2 V2 (w) ](4 ‘ $in *lf 


Comparisons have been made between the approximation (eq. 19) and the exact 
solution (eq. 7) at 0^ * 90° and 0 2 ■ 30° for both phase and absolute 
value of R. The comparison Is shown In table 1 *3r N ■ 16. The phase up 
to ai * 50,000 rad/sec varies within 2 # ; It starts oscillating beyond 

A 

u> ■ 20,000. The magnitude |Rj 2 (u>)| varies about 42 percent at w * 50,000, 
corresponding to a^k Q of 3.0, while at u, * 30,000 the error In magnitude 
Is only 5 percent corresponding to a^k Q of 1.82. Therefore, the approximate 
solution Is useful up to about ka 1 ~ 2. 

Equation 19 thus contains Information not only on the separation distance 
of the two rings (* 2 - ) but also on the radius of the rings a 1 and a 2> 

This equation permits the determination of the distance between rings in terms 
of the compactness behavior of the phase. Since the distance between U 2 - ) 

Is usually larger than either the radius a^ or a 2 , the first case of non- 
compactness as the frequency increases occurs at k Q (£ 2 - ) cos 0 2 * it, 

when A/2 * ; since k Q ■ 2 tt/A. Equation 19 can be written as 


(!) 

a*. Exact 

«*, Approximation 

ISBBI 

Approxlmat ion 

a, k 
1 o 

a 

‘VV k o 

Ba 

12.94 

12.94 

21,993 

21,993 

.06 

.097 

3.6 


39.11 

19.11 

80,706 

80,705 

.18 

. 29 

1.0 


65,60 

65.60 

89,513 

89,507 

.30 

. 48 

1.79 

7,000 

90.90 

90.89 

116,300 

116,280 

.43 

. 68 

2.51 

10,000 

125.20 

125.15 

123,622 

123,519 

. 6 

. 97 

3.6 

15,000 

137.084 

136.56 

48,901 

49,353 

.913 

D 

5.38 

20,000 

123.47 

124.37 

63,450 

63.695 

1.2 

H 

7.18 

30,000 

207.42 

206.83 

14,248 

13,584 

1.82 

2.91 

10.0 

40,000 

284.56 

286.52 

488 

441 

2.4 

3.9 

14.36 

50,000 

367.49 

367.89 

5 

3.5 

i 

3.0 

4.8 

17.9 


Table 1.- Comparison between exact solut ion and approxlmat ion In the 

limit of k a. < 1, at 0 ■ 90°, t' ■ 30*. 
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A(u>)e 


- IT < 0A £ V 


where 


R,W ■ r{ 0) - j[k 0 <i,M] J r, (2) U) 


and 


R 2 (u>) 


R^ (id) 


R< 0) - J^U)] 2 S< 2 >U) 


Ro(u>) * A(w)e 


lB(u) 


(20) 


The phase 0 Is not equal to zero when Q 1? and Q 1 0 are not even function 
of x. This Is a consequence of the non-zero convection velocity of the flow. 

A A 

For real R^w) and R^w) 

A 

0 Ri (iii) > R 2 (w) 

it R^ (u) < R 2 (w) 


( 21 ) 




The phase of the cross spectral density a Is 


° '('* ' 'i) 


tr ♦ 6 


It now becomes 


-fe^) 


ti Rj(w) > 


a Is undefined when Rj(a>) * (to) 


a ■ it ♦ 


(tfh^y V“> * «?<*>> 


Computational results of the phase a In terms of k^tp " cos 0 2* are 


shown In figure 9 by equation 14. The phase at all angles 0* Intercept at 

k Q ( l 2 * cos n a * {t — ““jp j w s * nce > Rjtw)* Th1s 


result thus provides additional Information toward unraveling the details of 
the source distribution from far-fleld measurements . 

3. Ring radius k a. > 1.- In order to simplify the hlqher frequency, 

m 01 

expression for large k Q a^ we use the following approximation for large N 


N -Ik a, cos 

! i 01 j 

J-i 


i r v -1k„a, cos ♦ 

H • d * 


n -k^ cos ♦ 




Since we have N ■ 16 or larqer, one can emoloy this approximation and 
equation 14 reduces to 


25 


a *1k COS I r* a a 

«1 lf2 (u.) • Ne 01 * [<),(«) ♦ 2Q 1 j U)cos(B 1 A)Jj o (A 1 - Bj ) 


♦ Ne 


1k l, cos 0- 

0 2 2 


[Qj ,2^' W ) * ? Ql t 2^' u, ^ cos ^ A 2 A ^o^ A 2 " 
^Q 2 (u») ♦ 2Q 2 (u))cos(B 2 A)] J q (A 2 - Ba) 


[q^w) + 2Qi 2(lj)cos( B^A ( ^1 ” B 2 ) 


(22) 


When k a, * 1 the power series exDanslon for J can be used to recover 
01 0 

the results of the previous subsection. 

For high frequencies, 1.e.» » 1 , we note from physical reasoning 

that 

l 2 - » a, 


* 

R 1 reduces to 


R^ ^ C^) - Ne 


‘Vl 


♦ Ne 


-V 2 


cos 02 


cos 0 2 




♦ 2Q 1 (u>)cos(B 1 A)] J Q (A 1 - B^ ) 

♦ 2Q22( w ) C0S ( b i a )] j 0 ( a 2 - V 
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Since the acoustic energy of high frequency decays with distance downstream 
from the jet exit, 

A A 

Qj(w) » Q 2 (w) > $ 2 , 2 ^ 


Therefore we have 

* -Ik.t. cos 0«»« a , 

R, 2 (u) . Ne 01 1 Q^w) ♦ 2Q, j (u)cos(B^ A)l J Q (A^ - B 1 ) (23) 

The magnitude of the cross spectral density function vanishes when 
», - B, ■ k 0 .,(1 - sin e 2 ) ■ v m 

where v Is the m-th root of J ( z ) 
m .0 

A A 

A plot of the normalized absolute value | 2 (w)|/R 1 ^u)! Is shown In 
figure 10. Note that the first two zeros occur at k Q a 1 ■ 2.4 and 5.52 when 
the magnitude of the f*r-f1e1d pressure reaches a minimum. At lower value of 
k Q a, the undulation of the coherence function Is due to cancellation 
resulting from the matching of the acoustic wavelength with distance between 
rings. This effect occurs at k Q a < 1 , a problem discussed previously. 

In order to confirm this cancellation effect the numerical results have 
been compared with results based on the first four zeros of the Bessel function 
of order zero, as obtained f rom taMes. The results of this comparison is 
presented In t.'.ble 2. 
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k o*l 

k o a l 

b> 

from table 

computed 

rad/sec 

2.405 

2.403 

38,300 

5.520 

5.520 

87,977 

8.654 

8.653 

137,900 

11.792 

11.790 

187.900 


Table 2.- Comparison of k Q 0 | calculated from zeros of the Bessel function 
J 0 (k 0 a) * 0 with results computed fr„m the monopole model. 

The agreement Is excellent, Indicating that the coherence function or 
the absolute value can be utilized In obtaining Information on the source 
radius. The radius of the second ring Is generally not directly 
available from the absolute value of the far-field pressure, but to Improve 
the signal one can search for it in the plot of the coherence. 

There are also other ways to determine the radius a ? . One Is by 
assuming a spreading angle between the two rings based on physical reasoning. 
Thus, from the far-field correlation the location and the diameters c p the 
rings, is well as the type and compactness of the sources, can be extracted. 
The results are given In table 3. In order to better assess the validity 
of the approximate solution, additional comparisons of phase and magnitude 
values were made for a large range of u>. 

CONCLUSION 

This paper contains the solution of a direct and an inverse problem 
on the radiation of sound from randomly fluctuating sources on two rings. 

This work concentrated on obtaining the properties of the signals originating 


u 

rad/ sec 

a 

exact 

50,000 

-116.017 

60,000 

- 30.875 

70,000 

54.004 

80,000 

138.863 

90,000 

43.721 

100,000 

128.579 

li 0,000 

-146.563 

i 20 , OOu 

- 61.705 

130,000 

23.153 

140,000 

- 71.990 

150,000 

12.868 

160,000 

97.726 

170.000 

-177.416 


a 

approx. 

i 

exact 

-115.989 

3.9003 

- 30.875 

2.6190 

54.004 

3.8417 

138.863 

1 .12398 

43.721 

5.6892 

128.579 

2.0669 

-146.563 

6.9836 

- 61.705 

5.3225 

23.153 

8.2173 

- 71.990 

2.2542 

12.868 

4.2972 

97.727 

7.60 F 7 

-177.416 

2.9407 


28 

|R 12 M 

approx. 

4.0821 

2.6448 

3.8017 

1.0762 

6.6684 

2.1127 

7.0143 

5.2844 

7.9814 

2.5394 

4.3947 

7.7301 

3.0221 


Table 3.- Comparison between exact and approximate solution: 
k Q a^ > 1 , at $2 * 0° , N * 16. 
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from a random medium with known statistical properties. From the Inverse 
radiation approach which utilizes the memory of the received signals 
associated with the distortion, several questions regarding physical Inter- 
pretation of the signals In relationship to Its origin have been clarified. 

In particular, the Interpretation of data related to phase and coherence 
functions permits the establishing of: (1) the center of mass of the singu- 

larities, (2) the geometrical shape and size, and (3) the ability to dis- 
tinguish one type of singularity from another and from it determine the 
basic directivity of the different components. In addition, from the signals 

one can (4) deduce the distribution of the singularities either when 
compact or noncompact with respect to the acoustic wavelength. 

The location of the center of mass Is an Important factor to establish 
In an experiment, otherwise error can easily result, inadvertently, especially 
at those frequencies in which the far-field distance becomes geometrically 
limited. 

The results presented are mainly related to the meridian plane. Addi- 
tional information comes from the simultaneous observations of both azimuthal 
and meridian planes with ^( T ) taken as an odd function of the time delay 
t. This Information together with the partial recovery of the functions Q's 
will not be discussed at this time. 
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